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Abstract 

Here we follow the basic analysis that is common for real and 
complex variables and find how it can be applied to a quaternionic 
variable. Non-commutativity of the quaternion algebra poses obstacles 
for the usual manipulations; but we show how many of those obstacles 
can be overcome. After a tiny bit of linear algebra we look at the 
beginnings of differential calculus. The surprising result is that the 
first order term in the expansion of F{x + 5) is a compact formula 
involving both F'{x) and [F{x) — F{x*)]/{x — x*). 
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1 Introduction 



We are very familiar with functions of a real or complex variable x which we 
can expand, in the mode of differential calculus, as 

F{x + 5) = F{x) + F'{x)5 + + . . . . (1.1) 

But what if we consider a quaternionic variable 

X = Xq + ixi + ix2 + fcxs, (1.2) 

where those quaternions i, j, k do not commute with one another. The small 
quantity 5 will also involve all those quaternions. How then can we expect 
anything as neatly packaged as Equation fll.ll) ? 

This is a long-standing challenge to mathematicians; and here we believe 
we have something new to offer. Previous work has attempted to extend the 
usual concepts of the derivative, see, for example, references [1], [2] and 
[3]. We shall, instead, focus only on the differential dF, as suggested by Eq. 
(II. ip . The method will be to start by examining the exponential function e^, 
in Section 3, and then use this result to build more general functions F{x), 
in Section 4. We succeed in finding a general formula for the first order term 
that is surprisingly compact, as shown in Eq. (14.31) : but the price we pay is 
that it is no longer local, involving [F{x) — -F(x*)]/ (x — x*) along with F'{x). 

In Section 5 we describe a geometric view of this situation; in Section 6 
we find the Leibnitz rule to be valid. Section 7 presents the second order 
term in the generalized expansion (11.11) and Section 8 contains a general 
discussion of these results. In Appendix B we show how this line of analysis 
can be extended to the non-commutative variables based on the Lie algebra 
SU(2,C). 

First, however, we make a brief visit to quaternionic linear algebra. 

2 Linear Algebra 

Suppose we have the linear equation 

ax + xh = c (2-1) 

where a, 6, c are given quaternions and we want to solve for the unknown 
quaternion x. If we were dealing with real or complex quantities, then all the 
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factors in Equation fl2.ip would commute and we would write x = {a + b)~^ c. 
But that does not work for quaternionic numbers and variables. 
In fact, the solution to Equation fl2.ip may be written as 

X = d~\ac + cb*), d = + a{b + b*) + b*b (2.2) 

where one recognizes that [a, d] = ad — da = 0. The symbol * means the 
usual complex conjugation, which changes the signs of all the imaginaries. I 
would imagine that this solution is not new but I do not know where to find 
that out. 

There is an alternative way of writing the solution: 

X = {a*c + cb)h~\ (2.3) 
which I will leave for the reader to explore. 



3 Exponential Function 

For a general quaternionic variable x, we can define the exponential function 
in the usual way 

= lim (1 + (3.1) 
and this leads us to the expansion 

f' ds e-'^ S e'^ + 0(5% (3.2) 
Jo 

This formula is well known to some people (I believe it is often credited to 
Richard Feynman but I do not know a reference); a derivation of it is given 
in Appendix A. 

If we now separate x into real and imaginary parts as x = Xo+ru^, where 
(yU^y = —1, and, furthermore, use the expansion 

= e^°(cos r + Ux sin r), (3.3) 

then we can write the expansion as 

= + gx^^ ^ ^ ^ ^^^^ s]+cUxSux + 0{6^)] (3.4) 

where the real quantities a,b,c are given by 

1 , sin2r , , 1 cos2r — 1 1 , sin2r , 
a=-(H ), b= , c=-(-H ). 3.5) 

2 ^ 2r ^' 2 2r ' 2 ^ 2r ^ ^ ^ 



In the special case where commutes with 5, then this reduces to the 
famihar formula for the exponential of a complex variable. 

We can now go on to the logarithm function, defined as the inverse of the 
exponential: 

e^=y, x = lny, 6^^+^) = y + A, x + 6 = ln{y + A). (3.6) 
Using the previous results, we eventually arrive at the expansion 

ln{y + A) = lny + Ay-' A + B[u,, y-^A] + Cu, y-'A + 0{A^), 

1 r 1 

A = -(rcotr + 1), B = C = -{r cotr - 1). (3.7) 

4 General Function F(x) 

For a general function F{x) of a quaternionic variable x = xq -^rux-, we start 
by assuming a representation as a Laplace transform: 

F[x) = Jdp f{p) eP"" (4.1) 

where p is a real variable. We then calculate an expansion for F{x + 6) using 
the formula derived above for the exponential function; and this leads to the 
following general formula 

F{x + 6) = F{x) + F'{x) 5 + ^ {-F{x) + F{x*)) [u,., 5\ + 

-^F\x) [u,,{u,,5]] + 0{5^), (4.2) 

where F\x) is the usual derivative of the function F{x) calculated as if x 
were a real variable. It is surprising how simple this formula appears. 
Another form of this formula (14.21) is 

F{x + 5)- F{x) = F\x) 6i + {F{x) - F{x*)) {x - x*)-' 62 + 0{6^), (4.3) 

61 = - {6 - 6 u^), 62 = ^(6 + uju^), (4.4) 

which will be discussed further in the next Section. 

If one considers the function F{x) = x", then the above formula, (14. 2p . 
is correct, as may be shown by induction. Thus it is also true for any power 
series, F{x) = J2n Cn x"'. 

It is an interesting exercize to show that the special formula (13. 7p for the 
logarithm function is in agreement with the general formula (14. 2p . 
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5 Geometric View 



The variable x, as written in Eq. fll.2p . may be viewed as a point in a four- 
dimensional space composed of the real line, for xq, and a three-dimensional 
Euchdean space, for the vector x = (xi, 0:2, xs). In the decomposition x = 
xq + rux, we see r as the length of the vector x = rx; and the unit imaginary 
quaternion Ux is the dot product of the unit vector x with the vector k). 

Now we want to recognize that the quaternion 6, which is used to displace 
the variable x when we write F{x + 6), can be given a different decomposi- 
tion in that four- dimensional space: as defined in the previous section, Eq. 
f l4.4p . 6 = 61 + 62. We can recognize that 61 is in the two-dimensional space 
composed of the real line and the radial direction along the vector x; and 82 
is in the two-dimensional space that is orthogonal to the direction of x. One 
can say, geometrically, that is parallel to x (perhaps writing it as 6\\) and 
62 is perpendicular {S±). Algebraically, this is represented by the relations, 

61 X = X 61, 62 X = X* 62- (5.1) 

This decomposition of 5 is a local process. It varies from one place x 
to another, rather like the familiar unit vectors in polar coordinates of two- 
dimensional Euclidean space. 

The authors of reference [T] have taken a similar approach, introducing a 
local unit imaginary (which they call iota) that is the same as what we have 
defined as Ux- However, they limit their differentiations to displacements 
that are restricted to this two-dimensional space: they allow only what we 
call ^1 without any of 62- In that way they merely reproduce what is known 
about ordinary complex variables. 

6 Leibnitz' Rule 

Let us now define the first-order differential operator V, from Eq. (14.31) . as 

F{x + 6) = F{x) +V F{x) +0(6'^) (6.1) 

with 

V F{x) = F\x) 5i + {F{x) - F{x*)) {x - x*)'^ 62. (6.2) 
Now we consider V{F{x)G{x)). It starts off easy: 

V{FG) = {F{x) + VF){G{x) + VG) - FG = FVG + {VF)G] (6.3) 
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but the next step is more interesting, substitution fl6.2p into 06.31) : 

V{FG) = F[G'{x) 6i + {G{x) - G{x*)) {x - x*)'^62] + 

[F'{x) 6i + {F{x) - F{x*)) {x - x*)-^62]G. (6.4) 

The question is, How do we move the quaternions 6 to the right, past the 
function G in the second term of Eq. fl6.4l) ? The answer is given by flS.ip : 
and we finally get, 

V{FG) = {FG' + F' G)6i + [F{x)G{x) - F{x*)G{x*)]{x - x*)-^ 62, (6.5) 
which corroborates Leibnitz' rule for this differential operator. 

7 Second Order Terms 

Let's return to the exponential function (13.1 p and proceed with the expansion, 

g(x+5) =e^[i+ ds e-'''6e'''+ dt ' ds e-^'^''^''6e'''6e'^ + 0{6^)]. (7.1) 
JO Jo Jo 

The better approach is to combine the exponential function and the 
Laplace transform from the beginning. Writing F{x + 6) = F{x) + F^^^ + 
F^"^^ + . . ., we will first re-do the calculation of F^^^ to see how easily it goes 
with the decomposition S = 61 + 62 and the algebra (15. ip . 

= J dpf{p) p ds e(^-^)P"(5i + ^2) e^P" = (7.2) 

J dpfip) p e^^' ds[6, + e^^'^^*-^') ^2] (7.3) 

and this leads immediately to the result Eq. fl4.3p . 
Now we look at 

= / dpfip) e^V dt * e-(*+*)f^(5e*P^(5e'P^ (7.4) 
J Jo Jo 

Again, we decompose 6 and after a bit more work arrive at the result for the 
second order term, 

= 1f"(x) 5l + (Fix) - F{x*)) {x - xT^ (52^1 - M2) + 

F'{x) {x - x*)-^ 662 + F'{x*) {x* - x)-^ 6261. (7.5) 
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8 Discussion 



One may ask what restrictions there are on the functions F{x) considered 
above. At first, one would say that they should be real analytic functions; 
having terms like xax where a is a general quaternion would certainly cause 
trouble. One can extend this condition slightly by allowing F{x) (but not the 
function G{x) in Section 6) to be a real function with arbitrary quaternions 
multiplying from the left. That is, the power series form F = J2n Cn x"' 
could have arbitrary numbers c„. 

This bias to the left-hand side can be reversed if we change the original 
steps fl3.2l) . setting s ^ 1 — s, and (14.11) . putting f{p) on the right-hand side. 

It is noteworthy that our differential operators are no longer local: they 
involve F{x*) along with F{x). Nevertheless, it is surprising how simple and 
how general the results obtained here are. 

The Taylor series we have discussed above are expansions about the origin 
X = 0. In the usual complex analysis such power series may be about any 
fixed point x = Xf] but such a quaternion constant put in the middle of our 
expressions would appear to cause trouble. That can be rectified by defining 
the new variable y = x — x f = i/q + Vy Uy and then using this new imaginary 
Uy to separate the displacement 6 = 6i + 62. 

One may also ask if this general method may be applied to some other 
kind of non-commuting algebra beyond the quaternions. I believe that some- 
thing very similar can be done starting with a Clifford algebra. Another 
example is given in Appendix B. 
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Appendix A 

Here we give a derivation of the formula (13.21) for any non-commuting 
quantities x and 6. 

N—m—1 r 

s n+^r-^ii+^r+oc*^)}. (A.2) 
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In taking the limit — * oo, we convert the sum over m to an integral over 
s = ^ and this yields 

= gx ^ ^ ^ Q(^2)^ (^ 3^ 

Jo 

Appendix B 

Here we shall extend the general method used above for a quaternionic 
variable to something built on a more general Lie Algebra - specifically SU(2). 
Here is the Lie algebra: 

[Jl,J2] = J3, [J2,J3] = Ju [J3,Jl] = J2, (B.l) 

where the three J's are understood to be matrices over the complex numbers. 
In particular we shall use the relations 

e^^' Ji e-^-^' = Ji cose +J2 sine, e^-^' J2 e-^^' = J2 cose- Ji sine, (B.2) 

which follow from (IB. ID . 

The new variable x is to be constructed with four real parameters as 

X = Xq I + Xi Ji + X2 J2 + X3 J3 (B.3) 

and we want to expand F{x + 6) = F{x) + F^^'' + 0{6^), where 5 is a small 
quantity in that same space of matrices as x. Our first step is to define a local 
coordinate system at the given point x. By a suitable linear transformation 
(rotation) of the Lie algebra we make the coordinate x appear as 

X = Xq I + rJs (B.4) 

where we recognize that r"^ = x\ + x\ + x\. 

We can now separate the displacement 5 = 5\\+ 5±_ as follows. 

8\\ = 601 + 63 J3, 6± = 6iJi + 62 J2- (B.5) 

Now we are ready to study the first order term in the expansion, again 
using the representation of F{x) in terms of the exponential function. 

= Jdp f{p) p e^^ ds e-'P^ 6 e'^^ (B.6) 
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Since S\\ commutes with x, the first part of this is simply F'{x) S\\. For the 
part with 6± we use the formulas (lB.2p . where 6 is replaces by —spr. The 
integrals over s are trivial and we merely write sin{pr) and cos{pr) in terms 
of e^*^'' to get our final result. 

F(x + 6) = F(x) + F'(x)6ii + {F(x + ir) - Fix - ir)} — 5± + 

{F{x + ir)+F{x-ir)-2F{x)] ^[J^,5,_]+ 0{5^). (B.7) 

It should be noted that the 5-related factors in Eq. (]B.7p can be written 
in the following way: 

[J3,^±] = -[x,<5], (B.8) 
r 

5± = -^ [x,[a;,5]], (B.9) 

5\\=5-5^. (B.IO) 

This means that we do not have to carry out the "rotation" that gave us Eq. 
(IB.4P explicitely; the talk about choosing a local coordinate system is merely 
rhetorical. 

I expect that this method can be extended to other Lie algebras, with the 
quantity 5^ subdivided into distinct portions according to the roots of the 
particular algebra. The system of Eqs. (IB. 90 . (1B.10|) would be adapted to 



make those separations, using the known values of the roots; and those root 
values would also appear in the final generalization of Eq. (lB.7p . 
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